Polarized ^He as a probe for short range spin-dependent interactions 
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We have studied the relaxation of a spin-polarized gas in a magnetic field, in the presence of 
short-range spin-dependent interactions. As a main result we have established a link between the 
specific properties of the interaction and the dependence of the spin-relaxation rate on the magnitude 
of the holding magnetic field. This allows us to formulate a new, extremely sensitive method to 
study (pseudo-) magnetic properties at the submillimeter scale, which are difficult to access by 
other means. The method has been used as a probe for nucleon-nucleon axionlike P,T violating 
interactions which yields a two-order-of-magnitude improved constraint on the coupling strength 
(Ssffp) as a function of the force range (A): g^gpX^ < 3 x 10~^^ m^. 

PACS numbers: 14.80.Va, 67.30.ep 



Hyperpolarized '^He is currently applied to a wide va- 
riety of scientific and medical problems. They include 
magnetic resonance imaging, spin-polarized targets, sur- 
face science, probing of biological systems, and precision 
measurements in fundamental physics |ll-llC| . One of the 
unique properties of polarized "^Hc is the very long spin- 
relaxation time constant which can be of the order of 
thousands of hours [ll|, 12 1, making polarized "^He ex- 
tremely sensitive to any spin-dependent interaction. It 
is well known that the presence of a magnetic field gra- 
dient in a cell containing spin-polarized gas significantly 
affects the spin-relaxation. The origin of this relaxation 
mechanism is the loss of phase coherence of the atoms 
due to the fiuctuating magnetic field seen by the atoms 
as they diffuse throughout the cell. Spin relaxation of a 
gas in the presence of a uniform gradient has been ana- 
lyzed by Gates, Schaefer and Happer using a perturba- 



tion theory approach [13[ and by McGregor [IJ] within 
the framework of Redfield theory. In the present work 
we focus on spin-relaxation phenomena due to the field 
gradient decaying over distances much shorter than the 
cell size. This results in new expressions for the relax- 
ation rates which agree with the already known result 
[13I [1^ in the limit of high pressure and high magnetic 
field (adiabatic regime) and in the uniform gradient limit 
for very low pressure and magnetic field ( "motion narrow- 
ing" regime) . A broad transitional regime of spin motion 
has been discovered which carries the complete informa- 
tion on the strength of this short-range field as well as 
on its spatial dependence. This finding allows us to pro- 
pose a new very sensitive method to study short-range 
spin-dependent interactions from various origins, e.g. the 
magnetic field of a strongly diluted ferromagnetic sam- 
ple or the pseudomagnetic field of hypothetical axionlike 
forces [ISMlTI I . To test the power of the method we have 
performed an experiment with polarized '^He measuring 
the longitudinal relaxation rate as a function of applied 
magnetic field. The experiment results in a new, stronger 
constraint on the axionlike interactions. 



Suppose the magnetic field in a cell may be described 
by a homogeneous magnetic field Bg and weak inhomo- 
gcneous field b(r): B(r) Bo + b(r), with < b(r) >= 0. 
According to Slichter [18[, from the Redfield theory of 
spin-relaxation due to a randomly-fluctuating magnetic 
field, the spin-relaxation rates are given by 
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[Sxi^l) + Sy{u))] , 



r2 = ^= ^[5.H + 5,H + 25.(0)] 
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Here Ti is the time constant for the longitudinal relax- 
ation rate, T2 is the time constant for the transversal 
relaxation , 7 is the gyromagnetic ratio for the atoms of 
the gas (7 « 2.04 x 10^ s'^ G'^ for ^Hc), and w = 7B0 
is the Larmor frequency. The functions Sk=x,y,z(j-^) are 
the Fourier transform components of the magnetic field 
autocorrelation function i?fe(r) = {bk{t)bk{t + r)), where 
the ensemble average can be evaluated as 

Rkit-to) = JJ p{ro,to)p{r,t\ro,to)bk{r)bkiro)d^r(fro, 

(3) 

knowing the conditional probability density p(r, i|ro,io) 
for an atom sitting at time to at tq to be found at later 
time t at a position r. The factor p{rQ,to) is the single 
probability density. For times |t — to| much longer than 
the mean time between atomic collisions, the conditional 
density obeys the diffusion equation [19[, with the con- 
straint that the initial p(r, tokoj ^o) = ^(''—''0) and reflec- 
tion boundary conditions. Let us consider a rectangular 
cell of length L with square base of size R. Let x be the 
cell axis so that the square ends occur at a; = zLL/2, a 
homogeneous fleld Bq directed along z axis, and an in- 
homogeneous magnetic fleld b(a::) directed along x axis. 
For this geometry the problem becomes unidimcnsional 
with a known analytical solution for p [l4] . Substituting 
this solution in ([3|) and taking into account that in our 
case p(ro,io) = 1/i we get the following expression for 
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the autocorrelation: 
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and 
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/■^/^ dx 
bx,n= b{x)sm{{2n + l)nx/L)— (5) 

J-L/2 ^ 



where we have introduced the characteristic diffusion 
time constant = I? jD , D is the diffusion coefficient 
of the gas. Taking the Fourier transform gives 
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To proceed further we have to specify the field b{x). To 



take a concrete example, we consider the macroscopic 
pseudomagnetic field representing an axionlike interac- 
tion of polarized "^He with the nucleons in the cell walls 

b{x) fe, (e-(^/2+-)A - e-(^/2-)A) , (7) 
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Ng^gp (1 - e-"/^) 



(8) 



Here, x is the distancec from the wall, gs and gp are di- 
mcnsionless scalar and pscudoscalar couplings between 
the nucleon and the axionlike particle, A = — ^ is the 
force range, to„ is the nucleon mass, N is the nucleon 
number density and d is the thickness of the walls. Sub- 
stituting ([7]) in (5)- (6) we arrive at a general expression 
for the longitudinal relaxation on the one-dimensional 
case: 



Q _L/An2 / (1-0a(0a+2)) sin ^0i,/2 + (1-0a(0a-2)) sinh ^4>l/2 \ 

Ti = (76a)VA ; ^ cos 7^r72+cosh 7^172 (Q) 

V^a + <f>l)' V+iv/0l72 (sechA)^ (1 + 02^ A(^2_3),i^l^L)y 




10" 10"' 10"^ 10' 1 10 10^ 10' 

Holding magnetic fieid B [G] 



FIG. 1: Relaxation rate due to a short-range gradient mag- 
netic field (O versus tlie magnitude of the homogeneous mag- 
netic field calculated according to (|9]) for three different values 
for the range A. 



where 4>l = ^t^, 0a = ujt\ and t\ ~ }? jD. 

Our result © is illustrated in Fig. 1. We can dis- 
tinguish three regimes of relaxation: relaxation in a low 
magnetic field to <^ with Tiocconst, relaxation in 
a moderate field r^^ <^ cj <^ t^^ with Tic^uj"^^'^ , and 
relaxation in a high field u 3> t^^ with Fiocw^^. 

The magnetic field value corresponding to the tran- 
sition between the latter two regimes depends on the 
force range A: the smaller A is, the higher magnetic field 
is needed. For an extremely sharp correlation function 



i?fc(T), like a Dirac delta- function, the Fourier transform 
5fe(a;) is constant and the corresponding relaxation rate 
is independent of magnetic field (relaxation due to dipolc- 
dipole interaction in atomic collisions [2l| or due to col- 
lisions with walls [22|) 
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A>L, (10) 
A<L. (11) 



The expression (fTOl) agrees, with numerical factors of 
order one, with that [l3|, [ij] obtained for a spherical cell 
in a magnetic field with a uniform gradient. The expres- 
sion (|11[) shows that in the "motion narrowing" regime 
the relaxation is governed by the mean squared value of 
the inhomogcneous field, and by the time required for an 
atom to diffuse over the cell. In the limit of very low 
pressure, when the mean-free path IcoW is larger than the 
cell size L, the natural time scale is tl = L/v, v being 
the particle velocity. With this replacement our expres- 
sion pT|) agrees with obtained for ultracold neutron 
depolarization. 

In the opposite "adiabatic" limit w ^ ^ 3> ^ , our 
result ^ reduces to 
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(12) 



uj^XL Bq dx 

also in agreement with [l3, 14|. For a finite-size rectan- 
gular cell and short-range X R, L axionlike interac- 
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tions with the cell walls, all components of the pseudo- 
magnetic field bx{x), by{y), and bz{z) need to be taken 
into account. Since the Brownian motions of the three 
coordinate directions are independent, from ([T]) follows 
Fi = Ti^x + Ti.yj where both ri^;^ and Fi^^ are given 
by our expression Q for the one-dimensional geometry. 
The transverse relaxation rate follows from eqns. ([2]) and 
reads, for a rectangular (or cubic) cell 

r2-Fi(0)==2(76a)VA for A«i?,i. (13) 

We expect similar result for other geometries (cylindrical, 
or spherical) since ([T^ is independent of the cell size. We 
checked by Monte Carlo simulation that this is indeed the 
case at the 10% level of precision. 

We now apply our theoretical results to search for an 
exotic short-range axionlike interaction. Assuming for 
the time being no extra interactions, the experimental 
spin-relaxation is determined by contributions from three 
sources: Tiexp = ^dd + ^waii + T^, where F^^ is the 
dipolc-dipole relaxation due to atomic collisions, T^aii 
is due to the ^Hc spin relaxation on the walls of the cell 
and Fm is due to magnetic field inhomogencitics. Tdd and 
^waii are expected to be independent of magnetic field. 
If Bq is high, so that uj ^ t^^, the relaxation due to 
magnetic field inhomogeneities can be decomposed into 
two components as follows: 

|2\ I /l„ |2\\ 



D{{\g.f) + {\g.f))/Bl 



(14) 



where gi is the gradient due to inhomogeneities of the 
holding magnetic field Bq and ge is the gradient caused 
by the external magnetic environment. Since giOcBg-, the 
relaxation F^e due to the magnetic field gradient caused 
by the environment is the only term that depends on Bq. 
As all external magnetic sources are far away from the 
cell containing the polarized '^He, the relaxation is only 
affected by magnetic sources which are large compared 
to the cell size, and hence providing a nearly uniform 
gradient over the cell volume. Thus, we can expect that 
F,„e scales as Bq^ . Finally, we can write: 

+ D\g,\^B^'. (15) 

The relaxation due to short-range spin-dependent forces 
^ depends very differently on the magnetic field com- 
pared to the simple law ^5]). and can thus be separated 
from the other sources of relaxation. To take advantage 
of this feature we performed measurements of the longitu- 
dinal relaxation of '^He as a function of magnetic field us- 
ing a cylindrical alumino-silicate glass cell (GE180, 5 cm 
diameter, 10 cm long, 3 mm wall thickness) filled with po- 
larized '^He gas (75% initial polarization at 0.3 bars pres- 
sure) placed inside a self-screening "Magic box" magneto- 
static cavity More details of the experiment will be 
published elsewhere. The experimental results obtained 
are shown in Fig. 2 together with a fit of expression ([T5|): 

Tdd,wall,m^ = (9.71 ± 0.16) X lO'^ h'^, 

\ge\ = (2.70 ± 0.05) X 10^2 Q/j^^^ (j^g) 
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JgJ =(2.7±0.1)x10"^G/m 
I d.o.f. = 0.8 
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FIG. 2: Experimentally measured longitudinal relaxation 
rate of polarized '^He versus magnitude of holding magnetic 
field (points). Solid line shows a fit of the theoretical predic- 
tion (15) for no short-range spin-dependent forces. 



Now we consider the addition of an axion rate (jH]) to 
the normal expression (|15p . The fit involves now two 
additional free parameters A and ba, it yields a result 
compatible with zero for ba- Thus the experimental data 
shows no evidence for a new axionlike interaction. In 
order to set an upper limit on the strength of an axion- 
like interaction, the fit was performed again with various 
fixed values of A ranging from one micron to one cen- 
timeter. For each such value of A an upper limit on the 
product gs<7p was derived from the fit, which is shown 
as bold dashed line in Fig. [3l A previous attempt [25| 
to determine constraints on axionlike forces in the range 
10~^ m < A < 10~^ m using available '^He relaxation data 



111 . |12| appears to be flawed. That analysis was based 
on expression (12), which is only valid in the limit of 
high pressure and high magnetic field: u) 3> t^^ ^ ''"l ^' 
which corresponds to A > 10^** m for the experimental 
conditions 11, 12| . 

During the time that the present studies were per- 
formed, new experimental data on the very long time 
constant of the transverse relaxation in the "motion nar- 
rowing" regime became available from the shielded room 
BMSR-2 in Berlin [2^ . We have performed an analysis of 
these data within our theory. For the experimental con- 
ditions of [2^ the transversal relaxation time constant 
may be written as follows [l^, [3] '■ 
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^ = ^ + ^ {S7BI + ^Bl + 2VBI) (17) 

where the first term represents the field-independent lon- 
gitudinal relaxation due to interaction with walls and 
atomic collisions, while the second term represents the 
relaxation due to long-range gradients of the magnetic 
field. The first term ([T7| was measured in a high mag- 
netic field: Ti = 85±5 h. The second term was estimated 
to be = 370 ± 64 h from the measured values of the 
gradients in (|17l) . The authors [2^ compared the sum 
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formulated, supported by a first experiment: scanning 
the relaxation rate Fi on the holding magnetic field Bq. 
In general, this method samples the inhomogcneous field 
b{x) at different spatial scales. Information about b{x) 
can be extracted from comparison of the data Fi (w) with 
the expectations of the model. In the particular case of 
the simple parametrization one can extract the am- 
plitude ba and the range parameter A. Our first exper- 
iment could be dramatically improved if performed in a 
"zero field room" such as BMSR-2, where the external 
gradient of the magnetic field is 4 orders of magnitude 
lower than in our "Magic box" . 

We are grateful to V. Nesvizhesky for attracting our 
attention to the axion problem and to E. Kats and R. 
Whitney for valuable discussions. One of us (A. P.) ex- 
presses his gratitude to C. Fu and T. Gentile for fruitful 
discussion during a visit at NIST. 



FIG. 3: Constraints to the coupling constant product of ax- 
ionlike particles to nucleons gsQp as a function of the range 
A of the macroscopic interaction. Bold solid line: from long 
T2 [2^ (present work); bold dashed line: from scanning the 
magnetic field magnitude (present work); dotted lines: from 
UCN precession and depolarization [28l. [29l|: from UCN grav- 
itational levels [sOl; from mercury spin precession [sij . 



of these two terms T2pred = 69 ± 4 h: with the experi- 
mentally measured T2exp = 60.1 ± 0.1 h and concluded 
that there was agreement between the theory and experi- 
ment. We can interpret the difference between these two 
relaxation time constants in terms of the effect due to 
axionlikc interactions, via (8) and (13): 



F2a — 
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= (2.1±0.8) X 10"^ h^i (18) 



A similar result was recently presented in [27[ . Although 
the mean-free path between atomic collisions was as large 
as ~ 8 X 10^"'' m, the limit ([T9| still holds even for 
a new interaction in the micrometer range. Indeed, the 
diffusion regime approach, in particular Eq. (13), is valid 
as long as the mean-free path is much smaller than the 
cell size. This delicate point has been checked by Monte 
Carlo calculations based on Eq. ([2]). The new upper limit 
(fT9|) is illustrated in FiglHwith the bold solid line. 

Spin relaxation of ^He provides a sensitive method 
to study sub-millimeter magnetic properties or axionlike 
interactions. Analysing recent measurements [26| per- 
formed in the best magnetically shielded room yields an 
upper limit on axionlike interactions almost 2 orders of 
magnitude better than that obtained from the published 
data on UCN spin precession [2^. A new method is 
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